Enhancement of superconductivity by local inhomogeneity 
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We study the effect of inhomogeneity of the paring interaction or the background potential on 
the superconducting transition temperature, Tc. In the weak coupling BCS regime, we find that 
inhomogeneity which is incommensurate with the Fermi surface nesting vectors enhances Tc relative 
to its value for the uniform system. For a fixed modulation strength we find that the highest Tc 
is reached when the characteristic modulation length scale is of the order of the superconducting 
coherence length. 



Many strongly correlated superconductors, and in 
particular high-temperature superconducting (HTSC) 
cuprates, exhibit inhomogeneous electronic and/or struc- 
tural phases at the nanoscale |3,l3il3l- The coexistence of 
HTSC and inhomogeneity suggests that the underlying 
inhomogeneities could be at least partially responsible for 
the high value of the superconducting transition tempera- 
ture. Emery and one of us proposed that HTSC is related 
to frustrated electronic phase separation, commonly ex- 
pected in strongly correlated systems . These ideas for 
inhomogeneous superconductivity have been further de- 
veloped in the context of stripes ja, • It is important to 
distinguish these and related scenarios for superconduc- 
tivity creation or enhancement by inhomogeneity from 
the conventional weak-coupling coexistence of supercon- 
ductivity and various density waves 0, IE S ^ • In the 
latter case, the density wave order inevitably suppresses 
superconductivity due to the competition for the Fermi 
surface electrons. 

It is therefore important to understand the nature 
of the interplay between superconductivity and inhomo- 
geneities. A complete description of the interplay is 
clearly impossible. However, in the case in which the 
characteristic energy scale responsible for the formation 
of the inhomogeneity is much larger than the supercon- 
ducting energy scale (the gap A), and where the residual 
interactions are weak, a description based on BCS the- 
ory should be reliable. The purpose of this work is to 
study the effect of such imposed inhomogeneity on su- 
perconductivity within the BCS framework. The origin 
of the inhomogeneity could be either electronic, as in the 
frustrated phase separation scenario, or structural, that 
is caused by local lattice distortions or non-uniform car- 
rier concentration due to doping irregularities. We will 
assume that these structures do not cause Fermi sur- 
face nesting either due to the lack of periodicity (e.g. 
random doping profile) or due to the periodicity being 
incommensurate with the nested momentum transfers 
(e.g. frustrated phase separation, or stripes). Under 
these conditions we generically find that inhomogeneity 
enhances the global superconducting transition temper- 
ature, Tr. At the mean-field level, the maximum is 



achieved when the characteristic length-scale of the in- 
homogeneities, L, is large, in which case the transition 
temperature is that of the regions with the highest lo- 
cal Tc- Upon including the effects of phase fluctuations 
we find that Tc is maximized when L is comparable to 
the superconducting coherence length ^ ~ vp/Tc- The 
increase of the transition temperature occurs at the ex- 
pense of the superfluid density, which is reduced in in- 
homogeneous superconductors relative to their homoge- 
neous counterparts. 

Inhomogeneous pairing: Mean-Field treatment. As a 
first example we consider a Hubbard model with an in- 
homogeneous attractive potential U{r) > 0, 
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where ^k = Ck — M ^^id ncr(r) = ctcr(r)ccr(r) is the occupa- 
tion number of electrons of spin a at position r. Within 
this model, our goal is to understand whether for a fixed 
average pairing strength U{r), a uniform or non- uniform 
U{r) yields a higher transition temperature, Tc- In the 
weak coupling limit, we can derive the BCS condition for 
the onset of superconductivity from the Hamiltonian , 
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where Aq = Ekp f^(k)(Cq/2-k/2+pTCq/2-k/2-pi), U{k) 

is the Fourier transform of the pairing interaction, and 
K(p) is the pairing kernel. The pairing kernel depends on 
temperature T and the mean-field (MF) superconducting 
transition is defined by the temperature at which the 
integral equation has a non-trivial solution. The kernel 
can be calculated from the normal state electron Green 
functions 0| 
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where Nf is the density of states at the Fermi surface, Vf 
is the Fermi velocity, T is temperature, and ln7 « 0.577 
is Euler's constant. Here we also introduced an explicit 
high-energy cut-off for the attraction, lud- For T ^ Vfp 
this expression reduces to the well known homogeneous 
result, K ~ N f \n[2 jlud/ (ttT)]. 

The modulation of the pairing interaction leads to the 
mixing between Cooper pairs with different center-of- 
mass momenta. For simplicity we first assume a har- 
monic modulation {Q = 27r/L) of the pairing, U{r) — 
U + C/q cos(Q • r). In this case the integral equa- 
tion (O reduces to a system of linear equations A„ — 

UK„A„ + {UQ/2)[Kn-lAr,-l+Kn+lAn+l] = 7W„™A„, 

where A„ = A(nQ -I- go) and iC„ = if(nQ -I- qo). The 
"parent" momentum qo defines the minimal momentum 
of a Cooper pair in the connected family A(nQ -f go)- 

The paring instability occurs at the temperature Tc, 
such that the largest eigenvalue of matrix M is equal to 1. 
In the uniform case, Uq = 0, this condition is UK{0) = 1. 
We will now prove that the mean-field transition tem- 
perature is greater than in the uniform case, Uq ~ 0. 
Consider the go = family and without loss of gener- 
ality, take Uq > 0. Since all of the matrix elements of 
A4 are non-negative, by Perron's theorem the max- 
imal eigenvalue is a positive number that is larger than 
any diagonal matrix element, including UK{0). Thus, 
generically, the superconducting onset temperature Tc is 
increased whenever Uq ^ Q. 

This result can be understood from an analogy with 
a quantum mechanical particle in a tight-binding chain. 
Defining a new variable A„ = UKnAn, the BCS condi- 
tion takes a simple symmetric form, A„ = {l/UKn)An — 
(t/Q/2C7) (A„+i + A„_i). The "hopping" term delocal- 
izes the particle and thus reduces the "energy" below its 
minimal on-site value (I/UKq). Clearly, this leads to a 
relative increase of Tc- 

Large Q limit. In this limit, the quickly oscillating 
coupling is ineffective at mixing different modes, so that 
the off-diagonal terms in A4 are rapidly decaying with 
n. We are then justified in keeping only a small portion 
of the matrix surrounding the n = term. The lowest 
order correction to the homogeneous result is obtained by 
considering couplings between Ao and A±i. The largest 
eigenvalue in this case is 
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Given the separation of energy scales, Tc ^ VfQ <C ujd, 
we obtain by solving A,nax = Ij 



Te = — Wi?exp[-l/7V/(i7 + 77)], 
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where 77 = Uq'^ Ki/[2{1 — UKi)]. While ry is positive, 
and since Ki log[uj]j/vfQ] decreases with increasing 
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FIG. 1: Critical temperature for the inhomogeneous neg- 
ative U Hubbard model with coupling U{x) = U + 
Uqcos{Qx). The thick line denotes the mean-field re- 
sult, where Tc,a = (27/7r)LJo exp[— l/A^/17] and T^^h = 
(27/7r)a;D exp[-l/iV/([7 + |(7q|)]. The dashed line shows 
the critical temperature once phase fluctuations of the or- 
der parameter are included. For <^ 1, the super- 
conductivity is first established locally in regions where 
U{x) is large, but macroscopic phase coherence is achieved 
at a lower temperature, bounded from below by Tc ; = 
{2-f/n)ujDeyip[~l/Nf{U ~\Uq\)]. 

Q, so does rj. For vpQ > ujd, Cooper pairs can no longer 
scatter off the quickly oscillating coupling landscape, and 
we recover the critical temperature for the homogeneous 
case. 

Small Q limit. In the limit <^ 1, the 

global MF transition temperature is determined by 
the regions with the strongest pairing interaction, 
Tc w {2"f/TT)LUDexp[-l/NfiU + \Uq\)]. The deviations 
from this result due to finite Q and the effects of the 
phase fluctuations are discussed below. 

Electron density modulation. Before going further, let 
us in parallel consider the case of homogeneous coupling 
C/(r) = U, with inhomogeneity caused by a background 
potential variation. In the simplest case of the harmonic 
modulation, the additional contribution to the Hamilto- 
nian (Q) is 



Hp = C^icrCia- COS Q • 
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It is easy to see that for particle-hole symmetric density of 
states (DOS), the linear in p contributions to the BCS in- 
stability condition equations vanish identically. For small 
values of Q, the modulation acts as a slowly varying shift 
in the local chemical potential with the amplitude pro- 
portional to \p\. Thus, we only get a linear in p contribu- 
tion for asymmetric DOS, N{e) = iV/ + N'{e - ep)- We 
then find that BCS equations are identical to the case of 
inhomogeneous pairing interaction with the modulation 
strength 
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Ginzburg- Landau analysis (Q^ <C 1). We now consider 
the general case of slow variation of the pairing strength 
and/or background potential. The Ginzburg-Landau free 
energy functional in the presence of inhomogeneity is 



F= -J drdr'K{r-r')A{r)A{r')+ J dr 
+a J drp(r)A(r)2 + ^ J drA(r)^ 
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Here we assumed that the order parameter remains real 
even in the presence of inhomogeneity. We include both 
the coupling of the superconducting order parameter to a 
density wave, as well as the inhomogeneity of the pairing 
interaction. For small amplitude modulation of the pair- 
ing interaction, U{r) = U + 6U{r) with |<^t/(r)| < U, the 
two mechanisms are formally equivalent. For particle- 
hole asymmetric DOS, from the above considerations, 
a = —UN'/Nf. For simplicity, we only consider the 
inhomogeneous U{r) case here. The pair susceptibility 
kernel is given by Eq. (PJ. In the long wave length limit. 



K{r - r') = S{r - r')Nf 
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where ^ = vp/T. Computing the variation of Eq. ^ 
with respect to the order parameter, we find the equation 
for a stationary solution A(r), 



-e2v2A(r) + 5(r)A(r) + -^A(r)3 = 0, (9) 
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As a first step, we determine the inhomogeneous mean- 
field (MF) transition temperature for the pairing inter- 
action U{y) = U + Uq cos(Q • r), with QS_ < 1. Close 
to the MF transition, the cubic terms in Eq. can 
be neglected. Expanding in Uq/U, and transforming 
Eq. © to Fourier space, we obtain a system of equa- 
tions connecting A(k) and A(k ± Q). For small Q, 
A(k ± Q) « A(k), and after expanding up to the sec- 
ond order in Q, we obtain 

- 5™a.A(k) = (ek)2A(k) - ^A'Q^dlA{k). (10) 

Here gmax denotes g{r) evaluated at U{r) — U + \Uq\ 
and A = \/UQ/{NfU"^) (note that there is no explicit 
constraint on the value of parameter A since it is a ra- 
tio of two small numbers). The MF transition temper- 
ature is determined by the smallest eigenvalue of the 
differential operator on the rhs. This eigenvalue corre- 
sponds to the "ground state energy" of a harmonic oscil- 
lator, ^QA/^/2. The corresponding transition tempera- 
ture T*^^ = Tma!^ exp {-^QA/V2), is only slightly less 
than the transition temperature, T^^^ for a system with 
a homogeneous pairing interaction Umax = U + \Uq\. 



More importantly, it is easy to see that in the limit 
of small Q£,, the order parameter is exponentially sup- 
pressed in the region of smaller pairing interaction rela- 
tive to its value at the peak, 



A„iin - A,naxexp {-A i^Q)-^} , 
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This, in turn, implies that phase fluctuations, which we 
discuss below, can reduce the global phase coherence 
temperature signiflcantly below T*^^. 

The expression of Eq. (|ll|l is only valid at the MF 
transition temperature, where A is infinitesimal. To de- 
termine A below T*^^ we need to solve the non-linear 
Eq. (Q. In a d-dimensional superconductor, with ar- 
bitrary smooth variation of U(t), the boundary of the 
"classically forbidden" region, g{r) > 0, is a (d — 1)- 
dimensional surface. Hence, near the boundary the prob- 
lem is essentially one-dimensional, and in the g(r) > 
region we can apply the standard WKB approxima- 
tion to solve the linearized Eq. 0. The prefactor is 
fixed by matching the WKB solution to the intermedi- 
ate asymptotic at the boundary xq, which can be ob- 
tained by solving the full Eq. © in a linear potential 
g{x) — g'{xo){x — xq)- We then find that in the partic- 
ular case of harmonic modulation discussed above, the 
order parameter distance d away from the boundary is 
approximately 



A(d) ^ T^/AQ^cxp \-A^{d/0^^^ 
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This expression is obtained assuming g'{xQ) ~ A^Q, and 
therefore valid only for the temperatures sufficiently be- 
low T*^^ ^ r*lt So long as T » T^'i {T^'i is the uni- 
form Tc of a system with pairing strength U — \Uq\), the 
distance from the turning point to the minimum point of 
A is d ^ L. Notice, that this expression depends on tem- 
perature not only explicitly, but also through £, — Vf /T. 

Phase fluctuation effects (still with <^ 1). K con- 
sequence of the large spatial variations in the mean-field 
A(r) is that fiuctuation effects are severe where A(r) 
is small. Of these, the most important fiuctuations are 
thermal fluctuations in the phase of the order parame- 
ter, i.e. where A(r) ~ | AM_F(i")|e*^ where Am_f(i") is the 
solution of Eq. (0, and 9{y) is assumed to be a slowly 
varying function of r. The free energy cost of such phase 
fluctuations can then be readily computed 



dvj{Y){\/ef 
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where the local superfiuid stiffness is 
J{v) = N}^^\Amf(j)\^ . In general, the phase or- 
dering temperature estimated using this as the effective 
Hamiltonian is reduced from the mean-field transition 
temperature (i.e. the temperature at which J(r) 
vanishes), but by an amount that depends on dimen- 
sionality, and on the spatial arrangement of the regions 
of suppressed stiffness. 
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For concreteness, we consider the case of a two- 
dimensional superconductor. At finite temperature, 
no true long-range order is possible |l3l |. However, 
at T < T/fT: binding of topological excitations into 
vortex-antivortex pairs leads to a state with quasi- 
long range order, which has a non-zero superfluid 
stiffness [lj|. While for homogeneous BCS supercon- 
ductors in 2D, the difference between MF and the 
Kosterlitz-Thouless (KT) transition temperatures is tiny, 
(yMF _ Tkt)/T^^ - T^^/Tf (where Tp is the Fermi 
temperature), for inhomogeneous superconductors, the 
suppression of Tkt, is generally much larger. For a 
smooth random distribution of J(r), an estimate of Tkt 
can be made based on the effective superfluid density, 

Tkt [1/J(r)r' « ^ JrmnJ-max (14) 

This expression has a particularly transparent meaning 
for the unidirectional "striped- like" variation of U{y) that 
we treated explicitly when solving the mean-field equa- 
tions, above. There, Jmax corresponds to the stiffness 
along the stripes and Jmin ^ perpendicular to the stripes. 
The corresponding anisotropic XY model directly leads 
to the result Eq. H14|) . In this case, we find 

T 

Tkt-t^T^A 

min (Tkt)- 

^KT 

Together with Eq. ^ for l^min{TKT) - A{L,Tkt), 
this equation implicitly defines Tkt- With logarithmic 
accuracy we find that Tkt ^ min(r*^^, w/Q/A). 

In any case, baring certain artificial geometries, it is 
clear that for a long wave length modulation, ^ 1, 
the Kosterlitz-Thouless temperature Tkt is exponen- 
tially lower than the MF transition temperature; on the 
other hand, for modulations with ~ 1, the phase fluc- 
tuation region is very narrow and Tkt ~ T^^^ . In this 
regime, the mean-field superconducting temperature is 
still exponentially enhanced relative to its value in the 
uniform state with the same average paring interaction 
strength, U. For even faster modulation, ^ 1, the 
MF transition temperature drops since the pairings in- 
teraction modulation averages out on the length scale of 
f . This trend is presented qualitatively in Fig. ^ 

For a "dirty" superconductor with a mean free path 
shorter than the clean coherence length, l = vpT < ^, the 
effect of phase fluctuations can be estimated in the same 
way as in the clean limit, with minor changes (involving 
prefactors) but with the coherence length redefined as 

Summary. We studied the effect of nanoscale inhomo- 
geneity on the superconducting transition temperature, 
Tc- We considered two possible kinds of inhomogeneity: 
the modulations of the paring strength and of the back- 
ground potential. In the weak coupling BCS regime, we 
find that inhomogeneity which is incommensurate with 
the Fermi surface nesting vectors enhances Tc relative to 



its value for the uniform zero center-of-mass momentum 
pairing. For a fixed modulation depth we find that the 
highest Tf. is reached when the modulation wavelength is 
of the order of the superconducting coherence length. For 
shorter wavelengths, the superconductor cannot take ad- 
vantage of the locally favorable conditions, while for the 
longer wavelengths, the global superconductivity is sup- 
pressed due to the phase fluctuations on the weak links, 
where the amplitude of the order parameter is signifi- 
cantly reduced. Although explicitly derived for s-wave 
superconductors, similar results will also apply to un- 
conventional superconductors in the presence of smooth 
(on the 1/kf length scale) inhomogeneities. Clearly over- 
simplified, the presented picture bears resemblance to 
the high-temperature superconducting cuprates, where 
considerable experimental evidencels'l indicates that the 
maximum Tc occurs at a crossover between a regime 
where Tc is controlled by the pairing scale and where 
it is a phase ordering transition. 
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